Oscillatory Driving of Crystal Surfaces: a Route to Controlled Pattern Formation 
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We show that the oscillatory driving of crystal surfaces can 
induce pattern formation or smoothening. The driving force 
can be of quite different origin such as a pulsed laser beam, an 
electric field, or elasticity. Depending on driving conditions, 
step bunching and meandering, mound formation, or surface 
smoothening may be seen in presence of a kinetic asymmetry 
at the steps or kinks (the Ehrlich-Schwoebel effect). We em- 
ploy a step model to calculate the induced mass flux along 
misoriented surfaces, which accounts for surface dynamics 
and stability. Flux inversion is found when varying the driv- 
ing frequency. Slope selection and metastability result from 
the cancellation of the mass flux along special orientations. 
Kinetic Monte Carlo simulations illustrate these points. 



PACS numbers: 05.65. +b, 05.70.Ln, 68.65.+g 

Sculpting surfaces at nanometer scale is of major tech- 
nological interest. Besides lithography, spontaneous pat- 
tern formation during crystal growth has been proposed 
as a tool to create large scale nanostructured surfaces. 
In this letter we point out that oscillatory driving ap- 
pears as an alternative route for pattern formation. We 
see two basic advantages in this method: first patterning 
and growth are separated, so that morphology is not a 
function of the growth process. Second it offers better 
control of the structure. An in situ and real-time control 
of the pattern becomes possible, opening a wide range 
of new applications. 

On the side of fundamental physics, oscillatory forcing 
has been a long standing source of intriguing nonlinear 
phenomena, such as the inverted pendulum problem [Q , 
and the Faraday instability |^| in fluid mechanics. Oscil- 
latory forcing also induces pattern formation in granular 
media ||. Moreover, parametric forcing of an ensemble 
of oscillators is known to lead, for example, to motion 
of domain walls In this letter, we show that oscilla- 
tory forcing of crystal surfaces induces macroscopic mass 
fluxes, leading to pattern formation or smoothening. 

In the past 15 years, a large number of studies were de- 
voted to surface roughening during crystal growth. We 
show that all the instabilities identified during growth 
(namely mound formation, step meandering, and step 
bunching) appear under oscillatory driving. The physi- 
cal origin of this effect can be related to ratchets, in the 
sense that kinetic anisotropy of steps (Erhlich-Schwoebel 
(ES) effect) is used to produce a mass flux along misori- 



ented surfaces, as pointed out by Barabasi et al [gj, who 
showed that AC electromigration should lead to direc- 
tional smoothening of surfaces. 

We first derive the mass flux along a misoriented (vic- 
inal) surface in order to analyse the stability of the sur- 
face. The main results are illustrated by Kinetic Monte 
Carlo simulations. We also briefly mention some exper- 
imental situations where our analysis applies, such as 
pulsed laser-induced mound formation on metal-vacuum 
surfaces, ultra-sound driven pattern formation on thin 
films [H, and pattern formation at metal-electrolyte in- 
terfaces subject to an oscillating electrochemical poten- 
tial. 

We consider surfaces where adsorption, desortion, and 
defect creation (such as bulk vacancies) are not allowed. 
Mass transport then only occurs through surface diffu- 
sion. Since the mean height of the surface does not vary 
with time, we have 



d t (h) = -V.(j) 



(1) 



where the brakets indicate that we have averaged over 
the timescales of oscillations, and 9t = d/dt. The mass 
flux (j) is the key quantity for surface dynamics. Step 
properties are supposed to be isotropic and direct step 
interactions via elasticity, or electronic surface states, are 
neglected. Steps then only interact via mass exchange. 
The surface flux has two components (j) = (J)n -I- (G)s 
where n = Vh/\Vh\ and s.n = 0. From Eq.(|l]) it follows 
that a nominal surface is stable (unstable) if (J) > 
(< 0). On a vicinal surface of mean slope mo ^ 0, an 
uphill flux ( J) > also brings a destabilizing contribu- 
tion to step meandering (i.e. undulations perpendicular 
to the mean slope direction). Nevertheless, in this case, 
the (G) component also intervenes in the stability cri- 
terion. A vicinal surface of mean slope m = \Vh\/a, 
where a is the lattice spacing, is stable with respect to 
step bunching (i.e undulation along the mean slope) if 
9 m (J) > 0. Otherwise it is unstable. 

As a first qualitative example of a mass flux induced 
by oscillatory driving, let us consider a vicinal surface 
for which the temperature oscillates alternatively and 
abruptly between two values. In the high temperature 
regime, the equilibrium concentration on terraces is high, 
and the ES effect is repressed. In the low temperature 
regime, the equilibrium concentration is low, and the ES 
effect is strong (i.e. adatoms cannot attach going down- 
step). We start from a regular vicinal surface (Fig. |l| (a)) 
in the low temperature regime. We first switch to the 
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high temperature regime: atoms detach from the steps 
and go on terraces (Fig. [j] (b)). Steps slightly retract 
to the left. Switching back to low temperature, adatoms 
go back to the steps but only attach from ahead (Fig. 
0(c)). The step goes back to its initial position. A net 
uphill flux results from the last part of the cycle (This 
could be interpreted as a many-particle " ratchet effect" ) . 




(c) 



FIG. 1. Mechanism for an uphill flux, see text. 



In the following we calculate the mass flux (J) along 
a misoriented (vicinal) surface, where steps are regularly 
spaced (with a distance t = 1/m). Dynamics of a vici- 
nal surface are decribed by the ID step model of Burton 
Cabrera and Frank S, modified by Schwoebel B in or- 
der to account for step kinetics. On terraces between 
steps, the adatom concentration only evolves via diffu- 
sion: 

d t c = Dd xx c, (2) 

where D is the adatom diffusion constant and dt denotes 
a time derivative. At the steps, mass conservation reads: 



— — = Do x 



Dd x 



(3) 



where is the atomic area, z denotes the step posi- 
tion, and the index + and — indicates the low and high 
sides of the step respectively. Relation (^) is valid when 
flc± <C 1, i.e. when the adatom concentration is much 
smaller than that of the solid. In order to describe at- 
tachment and detachment kinetics at the steps, the in- 
coming diffusion flux is related to departure from equi- 
librium M: 



Dd x c± = ±v±(c± 



^eq 



(4) 



where, v± are kinetic attachment coefficients. We define 
the kinetic attachment lengths d± = D/i/±, that are 
small for fast kinetics and large for slow kinetics. 
Let us first consider a sinusoidal pertubation: 



D = Do + D\ cos(woi) 
d± = d 0± + d 1± cos(u t) , 

Ceq = CQeq + C\ eq COs(wo^) 



(5) 



where quantities with index 1 are small and not neces- 
sarily positive. The mean flux going through a step is 



(J) = - — (d x c + + d x c. 



(6) 



To zeroth order in the perturbation the solution of Eq. 
(0) with the boundary condition ^ is c = c% q yielding 
zero contribution to (</}. The oscillatory nature of the 
perturbation yields zero first-order contribution. From 
the second order solution of Eqs. (|^,^) we obtain the 
mean flux: 



di + (ch — I 



(J 2 ) = 
Ado-sh) — 



n 



2 l + m(d 0+ + d -) 
di_(ch — 1 + Ad +sh) 



(1 + d 0+ d -X 2 )sh + A (d - + d 0+ )ch 



(7) 



where ch = cosh(A/m) and sh = sinh(A/m), and A 2 = 
iuJo/D. The mean flux (J2) results from a combination 
of oscillations of the equilibrium concentration and step 
kinetics. The frequency and slope dependence of ( J2) is 
in general complicated. We do not wish to be exhaustive 
here, but rather to highlight some important features. 
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FIG. 2. Hopping frequencies in a discrete model 

In the high frequency limit, where lvq is smaller than 
phonon frequencies lu p ~ 10 12 s _1 , (in order to avoid ef- 
fects such as stochastic resonance with phonons), but 
larger than atom hopping frequencies T, one can ques- 
tion the validity of Eq. (Q), which is based on a step 
model (Eq.(Q-[|)), in the "hydrodynamic limit" (i.e. low 
frequency). Typically, the slowest rate is that of de- 
tachment from steps r f» 10 12 exp(— E /k B T). Tak- 
ing typical values E ~ 0.5eV, T ~ 300^, one obtains 
Tq ~ I0 3 s~ 1 . We use a discrete ID model to check 
the high frequency limit lu p loq T . In this model, 
atoms hop to nearest neighbors -without interacting, 
on a frozen periodic vicinal surface of interstep distance 
£ = La. Dynamics is described with help of hopping fre- 
quencies Td, Tq , and Tf for diffusion, detachment and 
attachment on both sides of the step site (see Fig. 0). 
Hopping frequencies can be related to the parameters of 
the step model via D = QT d , d± = a(Tf /T d - 1) @, 
and flc eq = r^/r^. These relations can be inverted 
to express frequencies as a function of the BCF model 
parameters. Using Eq.(j5[) and in the loq — > 00 limit, av- 
erage frequencies can be used for atom hops. One then 
finds that detailed balance is broken to second order in 



2 



the perturbation: a steady mass flux along the surface 
is present. Since we considered a "tracer" model on a 
frozen surface, an additional condition is necessary to 
determine this flux, because the amount of matter on 
terraces is not fixed yet. We choose the coverage at the 
step site to be 9q = 1 ( i.e. there is always an atom at 
the edge of a step). We then find the following steady 
flux: 



^steady 



D flci eq m 



2(1 + m(d 0+ + d -)) \a + d 0+ a + d - 



'i+ 



di- 



Taking the limit uj — » oo in Eq. (fjj) leads to the same 
result, except that the a's in the denominator are absent. 
This suggests that the step model is valid in the high 
frequency limit as long as step kinetics are not too fast 
i.e. g?o± > ci- 
In the low frequency limit luq — > 0, the flux vanishes: 
as {J2)oo ~ cJq. In this limit, a wide variety of slope 
dependances can be obtained. We shall first focus on the 
occurence of flux inversion as frequency varies. When 
attachment-detachment is fast, i.e. when mdo± <C 1, 
(and for uiq — > 0), an expansion of (Q) provides: 



<■*»> = 



rtci 



eq 



48A) m 2 



\(d 1+ -d^) 



(9) 



By comparison to Eq. (||), it is seen that if (d]_/do- — 
di+/dQ+)(di- — di+) < 0, there must exist a frequency 
u>* for which (J2) changes sign. Hence, surface stability 
can be changed by tuning the frequency. 

In the limit of slow step kinetics mdo± 3> 1 (and luq — » 
0), we find: 



Qci eq CO 2 

2D rrf 



di+ 
do+ 



do- 



d 2 d 2 



(d 0+ +d -) 3 



(10) 



Comparing the sign of this expression to that of Eq. (|9|) , 
we see that (J 2) might change for a special slope m*. 
Two cases are possible: (i) If di+ < di-, and di + do- > 
di_c?o+, then a nominal surface is unstable, but the flux 
changes sign for m = to*. In this situation, slope se- 
lection is expected [0. (ii) If di+ > di_, and di+do- < 
di_do+, the surface is stable with respect to small fluctu- 
ations, but mounds of slopes larger than a special slope 
to* are subject to an uphill mass flux. A surface of 
finite initial roughness may be destabilized although a 
flat one is linearly stable: it is metastable. An analogous 
metastable situation occurs for morphologically unstable 
steps in presence of kink Ehrlich-Schwoebel effect during 
growth 0. 

Since mass fluxes are second order, they strongly de- 
pend on the temporal "shape" of the perturbation. As 
an example consider the low frequency square pertuba- 
tion mentionned earlier (Cf. Figj^). Here cos(woi) is 
replaced by sign [cos (u>ot)] in Eq. (g), leading to the fol- 
lowing flux, calculated for finite perturbation amplitudes 
(i.e. d±± and c\ q are not small): 



(J) 
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di + — g?i„ + 2m(cii+do- — di_do+) 
[1 + m(d + di)][l + m(d - di)\ 



(11) 



where do = do+ + do- an d d\ = di+ + di_. The slope 
and frequency dependance of this flux are different from 
the sinusoidal case [Eq. ([?])], leading to new conditions 
for surface stability. Taking do+ = d\ + = 0, c\ eq = 
co eq and di_ = do_ (case A), which corresponds to the 
case depicted in Fig [j], we obtain an uphill flux, leading 
to mound formation and step meandering. If instead 
di_ = —do-, a downhill flux is found (case B). Nominal 
surfaces should then be smoothened and step bunching 
is expected on vicinal surfaces. 
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FIG. 3. SOS simulations on a 256 x 256 lattice, greyscale 
represents surface height, luo/2-k — 0.1 MCSPS -1 , Eb = 1., 
E lb = 0.3, E s = 2, and T = 0.4. For (a) and (b) E u = 1.9; 
for (a) and (d) E\ s — —1.9. Starting from a flat nominal 
surface (a) is otbained after 5x 10 5 MCSPS. Mound formation 
is seen and r.m.s. roughness is w m 8. Starting from the 
patterned surface (a), we obtain (c) after 2 x 10 5 MCSPS, 
with w w 0.25. In (b) and (d) we have plotted the height 
minus that of the intial regular vicinal surface. Steps are 
initally parallel to the y axis. Meandering (b) and bunching 
(d) of steps take the form of ripples along x or y. 

We have performed Kinetic Monte Carlo simulations 
in order to show that the above mentionned analysis al- 
lows one to predict the surface response to an oscilla- 
tory variation of hopping rates. We use a simple Solid 
on Solid (SOS) model on a square lattice, where hops 
are accepted with an Arrhenius law, the activation en- 
ergy being equal to the bond energy Eb times the num- 
ber of in-plane nearest neighbors. An extra energy bar- 
rier E s , accounting for the Ehrlich-Schwoebel effect, is 
experienced during interlayer hops. Taking the lattice 
constant as the unit length, and one Monte Carlo Step 
per site (MCSPS) as the unit time, we have D = 1/4, 
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d+ = (i.e. d 0+ = d - = 0), d- = exp[E s /T] - 1. £ b 
and E s oscillate in a square fashion, between Eqi, — Em, 
and Eqi, + En, (and similarily for E s ) leading to case 
A or B depending on wether their oscillations are in- 
phase, or out-of-phase. We use the following parameters 
E b = 1., E u = 0.3, E s = 2, E u = ±1.9, T = 0.4, and 
ujo/2tt = 0.1. The resulting patterning of the surface 
correspond to the expectations, as seen in Fig. [|. 

We also studied the coarsening of these structures on 
a nominal surface. We have performed simulations for 
switching frequencies 1/10, 1/30, and 1/100. In all cases, 
surface roughness behaves as t@ , with (3 « 0.44 ± 0.02. 
The typical lateral length-scale on nominal surfaces also 
follows a scaling law t a , with a = 0.09 ± 0.02. The up- 
down asymmetry of the mounds is smaller than during 
growth. This might be related to the non-zero mean step 
velocity during growth |l(| . As expected, the surface can 
also be smoothened (see Fig. 3(c)). 

Let us now turn to some applications of this phe- 
nomenon. A pulsed laser could induce the abrupt tem- 
perature change leading to the uphill flux presented in 
Fig. |l|. Conditions similar to that of Ref. jlj] would be 
needed. For our analysis to be valid, the surface should 
not melt, and no dislocation should be induced. 

In recent experiments, high amplitude surface ultra 
sound waves were applied to an Al thin film ||, and 
pattern formation was observed. It is not clear wether 
this results from the oscillatory driving or from Grin- 
feld instability as pointed out in Ref. Moreover, fur- 
ther experiments would be needed in order to determine 
wether dislocations are present. 

The oscillatory driving of the potential in an electro- 
chemical cell can provide morphological changes. On the 
basis of the surface-embedded-atom-method (SEAM) 
|l3f we have surveyed the dependence of the ES barrier 
and the adatom equilibrium concentration on Ag(lll) 
as a function of the departure of the electrochemical po- 
tential from the potential of zero charge (pzc) JlJ] . We 
find that increasing the potential from to +0.85 V (rel- 
ative to the pzc) increases the adatom formation energy 
from 0.85eV to 0.92eV (which implies a decrease of the 
adatom equilibrium concentration), but also increases 
the ES barrier from 0.22eV to 1.16eV. This cycle thus 
should lead to case A. We also find that as the potential 
is cycled from to negative voltages, c eq still increases, 
but the ES barrier goes back up. Therefore cycling be- 
tween 0V and negative voltages we should be in case 
B. Hence, it should be possible to create or smoothen 
mounds depending on the chosen cycle. 

Note that oscillatory driving could also be applied dur- 
ing growth via the techniques described above, in order 
to smoothen or pattern the surface. 

In conclusion, oscillatory driving is a promising tool 
for in situ and real time control of the surface morphol- 
ogy. We have shown that it allows one to create and 
smoothen patterns on nominal and vicinal surfaces. A 



wide variety of behavior among which step bunching and 
meandering, mound formation, slope selection, metasta- 
bility, and mass flux inversion are found. Similar phe- 
nomena should result from a Kink-Schwoebel effect, as a 
consequence of non-equilibrium line diffusion, as shown 
for growth in Ref. jll| . 

In order to be more quantitative and to predict the 
typical wavelength of the instabilities, the next step of 
this study will be to consider stabilizing effects coming 
from line tension, step interactions, and nucleation. As 
in growth, initial stages of the instability on nominal sur- 
faces, as well as the shape of mounds crucially depend 
on nucleation. Including this effect is the main challenge 
for a global understanding of pattern formation via os- 
cillatory driving. 
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